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^ . Abstract 

1 ^ ' The extraction problem of information about the location and shape of the 

■ cavity from a single set of the temperature and heat flux on the boundary of the 
conductor and finite time interval is a typical and important inverse problem. Its 
one space dimensional version is considered. It is shown that the enclosure method 
developed by the author for elliptic equations yields the extraction formula of a 

^ , quantity which can be interpreted as the travel time of a virtual signal with an 

I arbitrary fixed propagation speed that starts at the known boundary and the initial 

Tlj" ' time, reflects at another unknown boundary and returns to the original boundary. 
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^ ' problem for the heat equation. Let be a bounded domain of R"(?t, = 1, 2, 3) with smooth 

■ boundary. Let D be an open subset with Lipschitz boundary of such that D C ^ and 
Q\D is connected. Let T be an arbitrary positive number. Let u = u{x,t) he a.n arbitrary 
non constant solution of the heat conduction problem: 

ut = Au in {Q\D)x]0, T[, 
Qu 

— = inaDxlO, T[, 
0) = inil \D. 

Here u denotes the unit outward normal vector field on d{VL\D) . The boundary condition 
for dujdv on dD means that D is perfectly insulated. Note that this is the simplest model 
of the heat conduction in a conductive body having a hole or cavity D . 
The problem to be solved is 



1 Introduction 

In this paper we consider one space dimensional version of a typical and important inverse 
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Inverse Problem. Extract information about the location and shape of D from the 
single set of the observation data u{x,t) and du/ du{x,t) on 9Qx ]0, T[. Note that T is 

fixed. 

The problem comes from the thermal imaging and the solution method may have applica- 
tions to, for example, the detection of corrosion. There are extensive studies of uniqueness 
and stability issue of Inverse Problem. In particular, it is known that the observation data 
uniquely determine general D itself under a suitable condition on the heat flux on dVl. See 
Bryan-Caudill [2] , Canuto-Rosset-Vessella [3] and references therein for more information 
about the issue. 

In this paper we are especially interested in seeking an analytical formula for the 
purpose and start with one space dimensional version of the problem. In one space 
dimensional case there is a way of calculating the so-called response operator from a 
general single set of the observation data. Moreover in Avdonin-Belishev-Rozhkov [1] it 
is shown that from the response operator one can extract the spectral data. Thus in one 
space dimensional case the problem can be reduced to the inverse spectral problem which 
has been studied well. However, the reduction is based on the boundary controllability 
for the heat equation and not an easy way. To our best knowledge there is no attempt 
for finding the extraction formula of information about the location and shape of D from 
the single set of the observation data without reducing to other inverse problems. In this 
paper, we present a direct approach using the idea of the enclosure method for elliptic 
equations introduced by the author [4] . 

So what is the enclosure method? Consider a non constant solution of the eUiptic 
problem: 

A-u = in \ D, 

— = o^dD. 
ov 

In [4] we considered the problem of extracting information about the location and shape 
of D in two dimensions from the observation data that is single set of Cauchy data of u on 
dVt. Assuming that D is given by the inside of a polygon with an additional condition on 
the diameter, we established an extraction formula of the convex hull of D from the data. 
In the method a special exponential solution of the Laplace equation played the central 
role. The solution takes the form g-^«g™ ('^+«'^ ) where r(> 0) and s are parameters; 
both uj and uj^ are unit vectors and satisfy u ■ = 0. The solution divides the space 
into two half planes which have a line {x\x ■ u = s} as the common boundary. In one 
part {x\x ■ UJ > s} the solution is growing as r — > oo and in another part {x\x ■ uj < s} 
decaying. We define a function I^^^^± (r, s) of the independent variable r with parameter 
s which is called the indicator function and can be calculated from the observation data: 

Jan ou ov 

In [4] it is clarified that the asymptotic behaviour of the indicator function as r — > oo 
depends on the position of half plane x - uj > s relative to D. Recall the support function 
hoioj) — supjjg jjX-uj. We say that uj is regular if the set {x\x-uj — hoioj)} n dD consists 
of only one point. Then we established: for regular uj /^^,^x(t, s)\s=hD{ixi) is truly algebraic 
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decaying as r — > oo. This means that: there exist positive constants A and such that 

^lim^T''|4,^x(T,s)|,=^^(^)| = A. (1.1) 
This fact is the core of the enclosure method. Since we have the trivial identity 

4,,x(T,i) = e--(*-^)4,,x(T,s), 

from (1.1) one could conclude that: if s > /i£)(a;), then the indicator function is decaying 
exponentially; if s = /i£)(c<j), then the indicator function is decaying truly algebraically; if 
s < /i£)(c<j), then the indicator function is growing exponentially. Moreover by taking the 
logarithm of both sides of (1.1), we obtained also the one line formula 

hm ^"g''-'-^^"'°^' ^/.^(a;). (1.2) 

T >00 J- 

By integration by parts (1.1) is equivalent to the statement: the integral 

-Thoiio) f 

IdD dv 



JdD OV 



is truly algebraic decaying as r — > oo. The key point for this fact is, roughly speaking, 
there does not exist a harmonic extension in a neighbourhood of the point xq with xo-ou — 

hoiuj). 

This is the essence of the idea of the enclosure method. However this is the case when 
the governing equation is elliptic [5, 6, 7]. How can one apply the enclosure method to 
non elliptic case? This is another motivation of the study and here we give an answer to 
this question. 

Now state our one space dimensional problem. Let a > 0. Let u — u{x, t) be an 
arbitrary solution of the problem: 

ut = u^x in]0, a[x]0, T[, 

u^{a,t) ^0 ioTt e]0,T[, (1.4) 

u{x, 0) = in]0, a[. 

Then the problem is: extract a from u{0,t) and Ux{0,t) for < t < T. 
Let c be an arbitrary positive number. Let 

v{x,t) ^e-'^^e""^ (1.5) 

where r satisfies r > and 



Z^-CT\l+ldl- — 



The function i> is a complex valued function and satisfies the backward heat equation 
Vt + Av — 0. Moreover e^^v has the special character 
• a s < cx + t, then lim,- >oo e'^^\v{x, t) \ — 
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• ii s > cx + t, then lim,- >oo e'^^\v{x, t)\ = oo. 

The half plane {{x,t) \ s > cx + 1} in the space time plays the same role as the half plane 
{x \ X ■ uj > s} for elliptic case. Changing c means changing normal vector of the line 

cx + t = s and corresponds to changing u. 

Definition 1.1. Given c > 0, s e R, define the indicator function Ic{t; s) by the formula 

7e(r; s) = e^' T {-v^{0, t)u{0, t) + m,,(0, t)v{0, t)) dt, t > c"" 
Jo 

where u satisfies (1.4) and v is the function given by (1.5). 
Our main result is the following extraction formula. 

Theorem 1.1. Assume that we know a positive number M such that M > 2a. Let T and 

T' satisfy < T' < T . Let c be an arbitrary positive number satisfying Mc < min {T, 2T'}. 
Assume that Ux{0,t) coincides with a polynomial oft on the interval ]0, T'[ one of whose 
coefficients are not 0. The formula 

hm i^g^£(liM ^ _2ca (1.6) 

T >00 J- 

and the following statements are true: 

if s < 2ca, then lim^ ,oo Ihi^', s)| = 0; 

if s > 2ca, then linv >oo \Icij', s)| = oo. 

This is an unexpected result. Under the condition on T integration by parts gives 

Ic{r; s)\s=ca e^'^ {vx{a,t)u{a,t) - U:c{a,t)v{a,t)) dt (1.7) 
Jo 

modulo exponentially decaying as r — > oo. Since the set ]a, oo[ x ]0, T[ and s — ca 
correspond to D and s = hoiuj), respectively and the right hand side of (1.7) corresponds 
to the integral (1.3) in the elliptic case as explained above, from the past experience we 
expected that this right hand side decays truly algebraically as r — > oo. If it is true, 
then we automatically obtain the formula 

y log|/c(r;0)| 

lim = —ca 

T *00 J- 

which corresponds to the formula (1.2). However, the fact is different from the expected 
value. See Figure 1 below. 
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Figure 1: (1.6) extracts the line cx -\-t — 2ca not cx + t — ca. 



One heuristic, however, technical explanation for this phenomenon is the following. One 
knows that u{x, t) can be extended to the domain ]a, 2a[ x ]0, T[ by the reflection x — 
2a — xaix — aasa, solution of the heat equation. Then integration by parts gives 

e'""' r {v^{a, t)u{a, t) - u^{a, t)v{a, t)) dt ~ e^^" {v^{2a, t)u{0, t) + u^{0, t)v{2a, t)) dt. 
Jo Jo 

However, the function e'^'^"'v{x, t) is exponentially decaying as r — > oo in a neighbourhood 
of the set {2a} x [0, T]. Therefore one concludes the exponential decaying of the function 
IciT; s)\s=ca- This argument suggests the existence of the phase function of IciT;s)\s=ca 
with a negative real part. Theorem 1.1 clarifies that it is essentially —rca. 

So the next question about (1.6) is: what is the value 2ca? One interpretation for 
this is: it gives us the travel time of a signal with propagation speed 1/c that starts at 
the boundary x = and the initial time t = 0, reflects by another boundary x = a and 
returns to x — 0. See Figure 2 below. 

This is a quite attractive interpretation. From this point of view the restriction T > Mc 
is quite reasonable since Mc gives an upper bound of the travel time 2ca. The idea behind 
this interpretation is the belief, solution of heat equation Uf = u,j.,j. can be considered as 
a suitable superposition of solutions of the wave equations Uu = -^u^x in an appropriate 
sense. This belief is coming from a well known fact for the solution of the initial value 
problem for the heat equation. Given initial temperature / = f{x), x G R" with compact 
support, let u = u{x,t) be a solution of the wave equation Uu = A-u for x G R", t > 
with initial values u{x,0) = f{x) and Ut{x,0) = 0. Then it is well known (e.g., see [8]) 
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Figure 2: (1.6) extracts the travel time 2ca of the signal with propagation speed 1/c 
started at the boundary x — 



that the function 

vix^t) =2 u{x, s)ds 

Jo v47rt 

yields a solution of the heat equation vt = Av for x G R", t > with initial values 
v{x, 0) = /(a;). Then the change of variable s = yields 

v{x,t)^^ / e-^'/^u{x,CVi)dC. 
Replacing i with t^, we obtain the quite impressive formula 

In one dimensional case, we have 

vix,t')^-^ / e-^'/'if{x + Ct) + f{x-Ct))dC. 

2 ^J'K Jo 

The point is: the function u{x,E,t) of the independent variables x, t satisfies the wave 
equation with propagation speed ^, uu = Am with initial values u{x, 0) = f{x) and 
Ut{x,0) ~ 0. This means that v{x,t'^) really contains a signal with an arbitrary prop- 
agation speed ^ — 1/c and can be obtained as a superposition of signals with several 
propagation speeds. If this is true also for a solution of (1.4), then it is reasonable to 
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expect that the observation data coming from the heat equation should contain some 
information coming from the wave equation with arbitrary propagation speed ^ = 1/c. 
This is another role of c. In fact, in the final section we will see that a similar formula 
to (1.6) is valid for the wave equation with propagation speed 1/c. This suggests that 
the indicator function for the heat equation is a mathematical instrument that picks up 
a signal coming from the corresponding wave equation with propagation speed 1/c. 

Note that, in Section 4 the reader will see that another, rather technical restriction 
T' > Mc/2 is redundant. However, following the discovery order, we keep the original 
statement since we think that the proof of Theorem 1.1 presented in Section 3 under the 
restriction and the homogeneous Neumann boundary condition is still interesting. The 
proof is a time dependent approach and is based on 

• a representation formula of the solution of (1.4) 

• the explicit form of the eigenvalues and eigenfunctions for the Laplacian with Neumann 

boundary condition 

• the asymptotic expansion of the special function 

oo -1 

^iH = E(-1)'t^— 3, Re«;>0 (1.8) 

fe=0 K +W 

which has been established by Olver [9]. Thus, one may think that the proof gives an 
application of the special function Si{w) to an inverse problem. 

In a forthcoming paper, we will study the corresponding problems for the heat equations 
with variable coefficients. 

It should be pointed out that having Theorem 1.1 does not mean the end of the inverse 
problem. Needless to say, the real data always contain an error. The next problem is to 
consider the case when the data are given by u{0,t) + Ei{t) and Ux{0,t) + E2{t) where 
the size of Ei and E2 are dominated by a positive number S. Here we mean by the size of 
an arbitrary square integrable function E — E{t) on ]0, T[ the norm of E over ]0, T[. 
We assume that 

\\{Ei, E2)\\ = ||^i||l2(0,T) + ||^2||l2(0,T) < S. 

Define the corresponding indicator function (at s = 0) by the formula 

7,(t; El, E2) = r i-v,iO, t){u{0, t) + E,{t)) + {u,{0, t) + E2{t))v{0, t)) dt, r > 
Jo 

where v is the function given by (1.5). 

We assume that 6 is known and sufficiently small. The following asymptotic formula 
(1.10) says that if one chooses a suitable r depending on 6, then the ratio 

logmT;E,,E2)\ 

T 

gives an approximate value of —2ca. 
Corollairy 1.2. Given 5 and a e]0, 1[ define 

T.{5) ^ U\og5\. (1.9) 
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sup I — h2ca|=0 — Yi — h (1-10) 



|lo 


g 


log 






log (5 





T/ien, as 5 — > the formula 

I log I 
sup I 

Needless to say, (1.9) is one of possible choices of r. In practice this choice of r will be a 
problem, however, it would be interested in doing a numerical testing of the formula (1.6) 
as done for (1.2) in [7]. This belongs to our future study. 

A brief outline of this paper is as follows. Theorem 1.1 is proved in Section 3. The 
proof is based on an asymptotic behaviour of an integral involving u on {a} x ]0, T'[ in 
the case when Ux{0,t) is given by a constant on ]0, T'[, as derived in Section 2. From the 
proof one may think that it is quite difficult to cover other boundary conditions at a; = a. 
However, in Section 4 we present a simpler approach that is based on the transform of u 
into w by the formula 

wixjT) = / u(x,t)e ^ dt 
Jo 

where 

1 



z"^ = T + i2cVi 



We see that the asymptotic behaviour of w{a, r) as r — > oo yields a corresponding 
formula for the cavity with the Robin boundary condition. We believe that this idea will 
work also for the multidimensional case. 

In the final section we give an application of this approach to the wave equation. 



2 A key lemma. Time dependent approach 

The integration by parts gives 

/c(r;s) = — e"/ u(a,t)vx(a,t)dt — e^^ uix,T)vix,T)dx. 
Jo Jo 

Since v{x,t) — e~^^*e^^, we obtain 

/^(r; s) = -ze^'e""' ^ u{a, t)e-'^*dt - e^e"^'^ T u{x, T)e''dx. (2.1) 
Jo Jo 



Rewrite this as 



Ur; s) = -^e^^e^"^ x e""^ / u(a, t)e-'^*dt 

Jo 



-ze^'e^"' X e-"' 



r uia, t)e-'^^dt - e'^'e-'^^ uix, T^e^'dx. 
Jt' Jo 



(2.2) 



Let s = 2ca. Since T and T' satisfy T > 2ca and T' > ca, respectively, we sec that 
the second and third terms of (2.2) have the estimates 0{Te~'^^-^ -c^)^ and 0{e~'^^-^ ^'^''"•^), 
respectively. These yield 



—2iacT\l 1 — 



Ic{r;s)\s=2ca^-ze V cH ^-az j u{a,t)e-'"'dt (2.3) 

Jo 
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modulo exponentially decaying as r — oo. 
Therefore it suffices to prove that 

rT 



f 2 

/ u(a,t)e~'' ^dt 
Jo 



is truly algebraic decaying as r — > oo. 

For the purpose we recall a known representation formula of the solution. 
Let fl be an arbitrary bounded domain of R". Let u satisfy 

ut = Au inOxjO, r[, 
u{x, 0) = in fl. 

The representation formula of u which is given below involves an associated elliptic prob- 
lem with the parameter t and eigenfunctions for the Laplacian with the Neumann condi- 
tion. This is taken from [2] . 

Let V — v{x,t) be the unique solution of the elliptic problem depending on t: 

1 f Ou 
Av = —— / —(x,t)dS(x) inO, 
\U\ Jan ou 

dv du , . 



/ 



vdx — 0. 



Let {Afcj^o {^k}'k=o be the all of eigenvalues and the corresponding complete 
orthogonal system of the Laplacian in Q with the Neumann boundary condition: 

A*fc + Afe*fc = infi, 



du 



~ ondQ, 



a, 



\-^k{x)\'^dx = 1 



and 

= Ao < Ai < A2 < • • • — > 00. 
Note that ^'^ with k > 1 satisfies 

/ 'ifk(x)dx = 
Jn 

since "^oix) = |ri|"^/^. 
Lemma 2.1 ([2]). The formula 

1 rK r du 



iff ou 

j^v{x,0)'^k{x)dx + j e~^''^^~^\j^vt{x,s)-^k{x)dx)ds^ ^fe(a;), 
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is valid. 

In the present case the domain fl is given by ]0, a[. v takes the form 

^'fe and Afe for A; > 1 are given by 

^k{x) = W-cos — X, Xk = ( — ) . 
V a a a 

A direct computation yields 



and 



/ v{x, 0)'i/k{x)dx = -u^{0, 0) J-Afe ^ 
Jo V Q 

f vt{x, s)'^k{x)dx = -Urot{0, s)J-Xl^. 

JO yd 



Therefore from Lemma 2.1 we obtain 

1 /■* 

u[x, 



i n 

,t) = v{x,t) / Ux{0,t)dt 

a Jo 

+\l-fl^k' (e-^'=*ii.(0,0)+ f\-'^>'^'-''>u,tiO,s)ds) ^k{x). 

V I 1 V ^0 / 



12 °° 

k=l 

Letting a; t a in this formula we obtain the expression of the boundary value of u at x — a 

1 /•* 



u{a,t) = -Ux{0,t) / u^{0,t)dt 

D a Jo 

-E(-l)%"Me"''^*«-(0'0)+ / e-'''^'-^^u,t{0:S)ds 
\ Jo 



2 ~ 

+ 



Now consider the simpler case 

u^{Q,t) = 1, Q<t<T' 

In this case we have 

a t 2 °° 



u^a, t) = ^-- + - J2{-l)'\,'e-'^\ 0<t<T'. (2.4) 



The following result is the key for the enclosure method and not trivial. 
Lemma 2.2. 

e-"^ / u(a,t)e-'^^dt 
Jo 

is algebraic decaying as r — > oo. 
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Proof. One sees that 
This gives 



Jo Xk + z^^ ' 



"^0 fe=i fc=i Xk{Xk- 



as T — )• oo modulo exponentially decaying. 
Moreover we see that 

and 



^2 



as r — > cxD modulo exponentially decaying. 
Using these and (2.4), we get 



/ u(a,t)e-'''dt - e-^'{— j + - V (-I)'^^tt ^} 



Here we cite the formula 

oo 1 2 

^ A;2 12 

«=i 

This yields 

Therefore the problem is the asymptotic behaviour of the function f{w) defined by 
the formula 

oo 1 

fe=l Xk + w 

We know that Olver (see p. 302 of [9]) has already studied the asymptotic behaviour of 
the function Si given by (1.8). Using the function one can write 



/(^) = -(^)^5,(^) + i^ 



' -i-rRii-) 



2^2 TT TT 
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where 

Then (2.5) becomes 

e— / u{a,t)e-^''dt - -^^R,{-^). (2.6) 

Here we make use of the formula (p. 304, (7.04) of [9]) for an arbitrary fixed S > 
and w — > oo in largu;! < — — S: 

^^^^^ - 2w- + r(l)«;V2 (, + l)-i/2 (2.7) 
where K_i/2 is the Macdonald function and has the asymptotic form 

K-i/2{z) = iY/^'e-^ (l + ^(;)) , ^ ^ oo in larg^l < ^ - 6. (2.8) 

See p. 250 and p. 238 of [9]. 

One imphcation of (2.7) and (2.8) is the asymptotic formula 

TTP-'^"' / 1 

R,(w)^- (l + 0(- 

w \ w 

This yields 



[ u{a,t)e-'''dt^-(l + 0{-)). (2.9) 



This is the desired conclusion. 

□ 



3 Proof of Theorem 1.1 and Corollary 1.2 

Consider the case when 

u^{0,t)^—,0<t<T'. 
ml 

We make use of the simple formulae: 

Ujn(x, t) = Ura-i(x, t)dt, m = 1, 2, • • • , < t < T' 
Jo 

where Um, fn — 0, 1, • • • are the solutions corresponding to the condition u^(0, i) 
t™/m!,0 < t < T. 

Using this formulae and the initial condition Mm(x, 0) — 0, we get 

Jo Z'^ Z'^ Jo 

z^ z^ Jo 
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and this thus yields the recurrence formula 



e-"' [ Um(a, t)e-'''dt = ^ / Um-i(a, t)e-'''dt + O fe-^^^'"^)) . 
Jo Jo ^ ^ 

This immediately yields 



-"M Um{a,t)e-'"'dt ^ —— uo{a,t)e-'"'dt. (3.1) 

Jo Jo 



Now wc arc ready to prove Theorem 1.1. 
By the assumption, one may assume that Uxi^it) has the form 

«.(0,t)= E^t™, 0<t<T'. 

Then the uniqueness of the direct problem gives 

I 



U{x,t) = E lraUm{x,t), < t < T' . 
m=0 

Prom (3.1) we obtain 

e-"^ r u{a,t)e-''^dt ^ V 7^6""^ Umia,t)e-''^dt 

Jo Z^n Jo 



m=0 



(3.2) 



Let j — min {m \ 7^ 7^ 0}. The number j satisfy ■jj ^ and 7^ = if m < j. 
A combination of (2.9) and (3.2) gives 

^lim^z'^+^e-"" £ u{a, t)e-'''dt = 27,- (^ 0) 
and from (2.3) we obtain the asymptotic formula 



2iacT\ 1 -— 

T >00 \ / I ij 

From this formula one knows that: there exist positive numbers Tq and A such that, for 
all T > To 

and 

^20+l)g2car|j^(^.0^|>^_ (3.4) 

Now all the conclusions in Theorem 1.1 follow from (3.3) and (3.4). In particular, we have 
log |/c(r; 0)1 ^ ^ 2(, + 1)^ + O f IV (3.5) 



T T \T 
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□ 

Finally we give a proof of Corollary 1.2. We have 

/e(T; E^, E2) - /c(t; 0) = T {-zEi{t) + E^it)) e-'"'dt. 

Jo 

This yields the estimate 

|/,(T;i?i,E2)-/c(r;0)| < (c^+-±=)s. (3.6) 



A combination of (3.4) and (3.6) gives 

\Ic{t;E,,E,)\ > |7e(r;0)| (l - ^e2-V(^+^)(cv/^+ (3.7) 

and 

\ISt-.E^.E2)\ < |7,(r;0)| (^1 + le'^W'^^+'\cV^ + -^)S^ • (3.8) 
Now choose r depending on 5 in such a way that 

2 2e<..^20-+l) + 1 )5 ^ (3.9) 



as 5 — > 0. Since we have assumed that T > 2ca, one possible choice of such r is to solve 
the equation 

The solution is just the r in Corollary 1.2. In this case, the estimation of the convergence 
rate of (3.9) is given by 

0(|log5|2(^'+i)+i/25i-'^). (3.10) 



Now from (3.5) and (3.7) to (3.10) we obtain the desired formula (1.10). 

□ 



4 A generalization. Stationary approach 

The aim of this section is to 

• present an alternative simpler method for the proof of Theorem 1.1 

• give a generalization of Theorem 1.1. 

First we describe the problem. Let a > 0. Let u = u{x,t) he an arbitrary solution of the 
problem: 

Ut = u^a: in]0, a[x]0, T[, 

Ua:{a, t) + pu{a, i) = fori e ]0, r[, (4.1) 

u{x,0) = in]0, a[ 
here p is an arbitrary fixed constant. 
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Then the problem is: assume that both of p and a are unknown, extract a from M(0,t) 
and u^{^, for < t < T. 

In this section the notion defined below plays the central role. 
Definition 4.1. Let c > 0. We say that a function / G ^^^(0, T) satisfies the condition 
c if there exist positive constant Ci, C2 and real numbers Tq{> c~^), Hi, such that, for 
all T > To 

Cit'^i < I / f(t)e-'''dt\ < Cst''^ 
Jo 



where 



2;^ = T + i2c^T^A 



It is easy to see that if f{t) is given by a polynomial of t on an interval ]0, T'[c]0, T[ that 
is not identically zero, then for all c > f{t) satisfies the condition c. Moreover, it is not 
difficult to sec that if f{t) is smooth on [0, T'[ and i = is not a zero of f{t) with infinite 

order, then f{t) satisfies the condition c. 

Definition 4.2. Given c > define the indicator function Ic{t) by the formula 

Ic{r)= r {-v^{0,t)u{0,t) + u^{0,t)v{0,t))dt, T>c-^ 
Jo 

where u satisfies (4.1) and v is the function given by (1.5). 

The following gives the answer to the problem mentioned above and generalizes The- 
orem 1.1. 

Theorem 4.1. Assume that we know a positive number M such that M > 2a. Let c be 
an arbitrary positive number satisfying Mc < T. Let the Ux{0,t) satisfy the condition c. 
Then the formula 

hm i^sH^ ^ _2,,^ (4.2) 

T ►OO 7- 

is valid. 

Proof. Introduce a new function w by the formula 

rT 



f 2 

w(x) = / u(x,t)e~^ ^dt, < X < a. 
Jo 



Note that, for simplicity of description we omitted indicating the dependence of w on r 
and c. This function satisfies 

w — z^w — '^u{x,T) in]0, a[, 

w'{a) + pw{a) — 
and the indicator function has the expression 

/^(r) = -^w(O) + w'(0). 
Then the integration by parts gives another expression 

/e(T) ^-{z + p)^/;(a)e"^ - e"^'^ T u{i, T)d^. 

Jo 
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This yields 

7,(r)e-'"^ = -{z + p)w{a)e-''' + O (e-^(^-'^«)) . (4.3) 
Now we write w{a) by using ^'(0). It is easy to see that w has the form 

w(x) = Ae^^ + Be-""' + f e"" u(^, T)e-^'d^ - e"^^ T u(^, T)e^'dA 

2z \ Jo Jo J 

where A and B are constants to be determined later. Then the boundary conditions at 
X — ^ and X — a give the system of equations 

zA- zB = w\^), 

{z + p)^e"^ -{z- p)Be-''' 



(1 + ^)e'^^ r ui^, T)e-^'di + (1 - ^)e-"^ /" T)e«^de) ■ 
z Jo z Jo J 



2 

Solving this system of equations, we obtain 

1 



z ((1 - ^)e- - (1 + ^)e-) 
X f(l- V(0) + ^{(1 + /\(e,r)e-^^de+(l--) r^^(e,7^)e«^rfe}) 



zf(l-^)e-"^-(l + ^)e" 
V 2; 2; 



X ( (1 + ^)^'(O) + ^{(1 + ^) r«(e,T)e-«^d^+ (1 - ^)e-2- P u{^,T)e^^dO 
^ z 2 z Jo z Jo J 

A direct computation yields 

2w'(0) 



i(;(a) 



z ((1 - ^)e-"^ - (1 + ^)e'^^) 



^ ^ fru(^,T)e^^d^+ ru(^,T)e-^^dA 



and thus this yields 



{z + p)w{a)e-''' - 2(1 + ^ + O(^))w'(0) (4.4) 

Z T'^ 



modulo exponentially decaying as r — > 00. Prom the assumption on Ux{^it\ (4.3) and 
(4.4) one concludes that, there exist positive constants C(, C'2 and real numbers Tq{> c~^), 
P\ and P2 such that for all r > Tq 
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This gives (4.2). 
□ 

Remcirk 4.1. If once a is known, then one can extract p by the formula: 

hm ^ ^ ^ ^—^ = -2p. 

5 Remark. Application to the wave equation 

Finally we give a comment on an application to the wave equation. Let a > and c > 0. 
Let u — u{x, t) be an arbitrary solution of the problem: 

utt = \u^x in]0, a[x]0, T[, 
-ii,(a,i)+pii(a,i) =Ofori e]0, r[, 

c (5.1) 
u{x^ 0) = in]0, a[, 

Mt(x, 0) = in]0, a\ 

here p is an arbitrary fixed constant. The quantity 1/c denotes the propagation speed 
of the signal governed by the equation. In contrast to the previous sections, one can not 
change c. 

Then the problem is: assume that both of p and a are unknown, extract a from u{0, t) 
and Mx(0,t) for < t < T. 

Introduce the function w by the formula 

fT 

w(x) — / u(x,t)e ^^dt, Q < X < a. 
Jo 

Again we omitted indicating the dependence of on r. Then w satisfies 

—w — T^w = e^^^luAx, T) + Tuix., T)) inlO, a[, 

-w'{a) + pw{a) = 0. 
c 

Now it is easy to see that 

2w'{0) 



w{a) 



cr(^(l + ^)e'=«^-(l-^)e-^"^) 



T ( (1 + ^)e^'^^ - (1 - ^)e- 

T T 



-car 



(5.2) 



X 



£{ut{^, T) + ™(e, T))e-^«^cie + /J(«t(e, T) + ™(e, T))e'^x) ■ 
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Definition 5.1. Define tfie indicator function I{t) by tlie formula 

^i^) = r f--^=^(0, t)u{Q. t) + -u^{0, t)v{0, t)) dt, r > 
Jo \ c c / 

wliere u satisfies (5.1) and v is tlie function given by 

v{x,t) = e-"(^^+*). 
Integration by parts gives tfie expression 

/(r)e2-- = (r - p)w{a)e^^^ - ce-^^^'"^^ /"(«,({, T) + Tu{i, T))e-^^dC. 

Jo 

Hereafter, using a similar arguments as done in Section 4 and (5.2), we obtain the following 

Theorem 5.1. Assume that we know a positive number M such that M > 2a. Let Ux{0, t) 
satisfy the condition: there exist positive constants Ci, C2 and real numbers tq{> 0), /xi, 
^2 such that, for all r > tq 

CiT^' < I r uJO,t)e'^'dt\ < C2T^\ (5.3) 
Jo 

Let T > Mc. Then the formula 

hm ^ _2ca, (5.4) 



T *00 



is valid. 

Needless to say, if once a is known, then one can extract p by the formula: 

cr(/(r)e2-''" + -w'{0)) 
hm — ^ = 2p. 

00 2w'(0) 

The quantity 2ca coincides with the travel time of a signal governed by the wave 
equation with propagation speed 1/c which starts at the boundary x = and initial time 
t = 0, reflects another boundary x = a and returns to x = 0. The restriction T > Mc is 
quite reasonable and does not against the well known fact: the wave equation has the finite 
propagation property. The condition (5.3) ensures that Ux{0, t) can not be identically zero 
in an interval ]0, T'[c ]0, T[. Therefore surely a signal occurs at the initial time. However, 
it should be emphasized that the formula (5.4) makes use of the averaged value of the 
measured data with an exponential weight over the observation time. This is a completely 
different idea from the well known approach in nondestructive evaluation by sound wave: 
monitoring of the first arrival time of the echo, one knows the travel time. 
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